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ABSTRACT - 

1 

t 

This memorandum deals with the capability of the Skylab 
CMGs to acquire and hold an earth-pointing (Z-LV) attitude. In 
this attitude experiments m the Earth Resources Experiment Pro- 
gram (EREP) are performed. The results of an MSFC study on this 
subject are reviewed and used as a point of departure. Included 
herein are: 


1. an operational procedure for obtaining the value of 
the total CMG angular momentum vector, H, at an appropriate 
point m the orbit, that minimizes the amount of TACS propellant 
consumed in acquiring and holding the Z-LV attitude; 

2. an explanation of the way in which this biasing of H 
can be achieved by existing ATMDC software; and 

3. some preliminary observations on a method for extending 
successive orbit capability in the Z-LV attitude with no concomitant 
increase m TACS propellant consumption. 

A comparison of the amount of TACS impulse consumed 
if H is appropriately biased with the impulse required by the 
MSFC approach is presented in Table 1. All of the results shown 
are for 120 deg experiment pass regions that are symmetrical with 
respect to orbital noon. 

TACS IMPULSE CONSUMPTION 
FOR BIASED AND UNBIASED H 



O 

o 

11 

o a 

8 = 30° 

8 = 50° 

UNBIASED H 

20 lb-sec 

28 lb-sec 

24 lb-sec 

BIASED H 

8 lb-sec 

8 lb-sec 

4 lb-sec 






TABLE 1 


JL 


* 



BA— 1 46 (8-64) 


BELLCOMM, INC. 


™- 70-1022-21 


DISTRIBUTION 


COMPLETE MEMORANDUM TO COMPLETE MEMORANDUM TO (CONTINUED) 


Martin-Marietta 

E. F. Bjoro /Washington 

Be 11 comm 

# 

A. P. Boy sen 
J,'. P . Downs 
, D. R. Hagner 
W. G. Heffron 
D. P. Ling 
J, Z. Menard 
P. F. Sennewald 
J. W. Timko 
R. L. Wagner 
M. P. Wilson 

Departments 2031, 2034 Supervision 
Division 101 Supervision 
Department 1024 File 
Division 102 
Central File 

P. R. Kurzhals/AMPD t 

MSC 

W. R. Cunningham/CB 
R. L. Frost/KS 
0. K. Garriott/CB 

F. C. Littleton/KM 
R. M. Machell/KW 
0. G. Smith/EW 

H. E. Whitacre/KM 

MSFC | 

W. B. Chubb/S &E-ASTR-SGD 

G. B. Har dy /PM- AA- E I 

H. F. Kennel/ S &E-ASTR-A 

G. F. McDonough/S &E- DIR 
E. F. Noel/S&E-ASTR-SI 
J. W. Sims/S&E— ASTN-PTA 
R. G. Smith/PM-SAT-MGR 
J. W. Thomas /PM- AA- El 

H. E. Wor ley /S &E-AE R0- DO 




Bellcomm. Inc. 



The TACS impulse capability is 61,000 lb-sec of 
which 30,259 lb-sec are allocated for performance of the 
total mission [2] . Included in this allocation is 100 lb- 
sec per pass for each of 45 Z-LV passes. The 30,741 pound- 
seconds of the total capability that is not allocated is a 
margin for contingencies such as operation with one failed 
CMG. In view of the results in Table 1 for unbiased H, 45 
Z-LV passes can be performed without biasing and stilT 
provide a good margin over what is now allotted. However, 
the results of Table 1 for biased H show that biasing would 
allow an increase m the number of Z-LV passes (without 
exceeding any margin requirements) should that be desired. 
In addition, since the TACS firings in [1] as well as this 
study occur during the experiment pass, a reduction in the 
number of firings is beneficial to experiment pointing 
stability. 


)l 
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TECHNICAL MEMORANDUM 


Introduction 


A recent MSFC study [1] explored the possibilities 
of biasing* the total CMG angular momentum vector, H, such 
that the maneuvers associated with the performance of the 
EREp could be achieved by the nested control law with a 
minimum expenditure of TACS propellant. The main conclusion 
was that, while some TACS propellant could be saved by such 
biasing, not enough could be saved to warrant any increases 
m ATMDC software complexity. The constraints imposed on 
the study were: 

1* | $ | <. 30 ° , and 

2. all experiment passes symmetrical about orbital 
noon. in addition, some preliminary results were given 
for 2 CMG operation and a statement was made that the 
spacecraft could recover from a Z— LV pass during the next 
CMG dump cycle with no TACS expenditure. 

The purpose of this memorandum is twofold- First, 
since [1] contains only results and conclusions, it is felt 
that a report that explains the procedure and the rationale 
behind it is necessary. Secondly, since the publication of 
[1] , certain changes in the planned operating procedure of 
the EREP have been either suggested [2] or made, e.g. the 
opening of the constraint on the solar incidence angle to 
|3| <. 50° [3]. The possibility of a more ambitious EREP also 
suggests a reconsideration of the need for a recipe for 
biasing H, especially if it could be accomplished with no 
increase in the software of the control system. Therefore, 
some extensions of the results of [1] are* also presented 
which, while not optimal, represent improvements in system 
performance with no increase in on-board software. 


*The aptness of this term will become apparent in the 
sequel. 
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II. A Z-LV Acquisition and Hold Procedure 

Figure 1 illustrates a Z-LV acquisition and hold 
procedure that has been synthesized from [ 1 ] and from the 
Program Definition Document (PDD) [4] . The goal of this 
procedure is to maneuver the spacecraft to an earth-pointing 
attitude in such a way that: 

1 . experiments may be performed during any 60° segment 
of the 120° region centered on orbital noon. (The variable 
Y 3 is employed in Figure 1 to permit consideration of 

unsymmetrical experiment pass regions) ; 

2 . the entire maneuver sequence can be executed by 
the nested control law rather than requiring TACS only 
control; 


3. the impact on the thermal and power systems is 
minimal . 

A complete description of this procedure requires some pre- 
liminary information concerning the momentum sampling strategy 
contained in the dump maneuver portion of the PDD.* Only those 
details that are necessary for a logical development are included 
in this section. Further details are in Appendix A. 


The orbital positions designated as ^ and 

(o) in Figure 1 are related through the characteristics of 
the sinusoidal component of the gravity- gradient torque, which 
is the most predominant disturbance on the Skylab spacecraft 
in the solar inertial attitude. This is most easily seen by 
considering the case when 3 = 0 ° and when the Y-geometric axis of 
the spacecraft is parallel to the normal to the orbital plane. 
For this case the gravity-gradient torque ,Ty, acting along 

the Y-geometric axis is as shown in Figure 2 where it is 
plotted versus orbital angular excursion, 0 . Using MSFC 
inertia data [5] , the analytical form of the curve in Figure 2 
is 


T ( 0 ) = -7.085 sin2{0 + a). 


where 0=0 is defined at orbital noon of the Z-LV pass , and 
the offset angle, a, is due to the fact that geometric, not 


* Because the primary spacecraft attitude for Sky lab is solar 
inertial, the momentum dump scheme is designed to work for that 
attitude. However, as shown in Appendix A, this scheme can be 
used as well to bias H for orbits that include Z-LV passes. 
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principal, axes are kept aligned with the solar inertial 
attitude by the control system.* One possible plot of the 
CMG momentum vector variation due to (e) is shown in 

Figure 3. The analytical expression corresponding to this 
plot is now derived. 


H y (t) = H y (t„) - j 7.085 sin2 ( w 0 t+a)di , 

to 

where the fact has been used that 0=ui o t, 0 ) o = the orbital 
angular rate, and t 0 is some initial time. The result of 
performing the above integration is 


H y (t) = H y (to) + 3167 [cos2 (tOot+a) - cos2 (w 0 t 0 +a) ] . (1) 

The plot in Figure 3 corresponds to the selection of 


H y (t 0 ) = 3167 cos2 (io 0 t 0 +ct) . 


(2) 


Equations (2) and (1) in combination with the variable changes 
e=u> 0 t and 6 0 =io 0 t 0 lead to (3) . 


H y (0) = 3167 cos2(0+a). 


(3) 


The adjustment of the level of the curve of H (6), e.g. by 
selecting H y (0 o ) for a given e Q is called biasing. 

Orbital positions ( 5 ) , (6) and ran now be 

defined. With respect to (3) , dump midnight, (jp , is 
located by (0+a) = +ir. The reason for specifically denoting 

is that all sampling points inthe ATMDC momentum dump 
strategy are measured relative to . Points Qp and 

(0) are located, respectively, by 0+a = -5tr/4 and 0+a = -3 tt/4. 


*The baselined Z-axis rotation, which essentially puts the 
X-principal axis in the orbital plane, is neglected. Its effect 
here is very small. 
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At these points is equal to its average {or bias) value 

(in this case, zero*) . Note that this is true only if no 
other torques except T^. in Figure 2 have acted on the Y 

axis of the spacecraft and if the spacecraft was in solar 
inertial for the entire daylight orbit proceeding Qf) . 

If the spacecraft were in solar inertial for the orbit 
preceedmg and if a bias momentum were built up during 

that orbit, tKen this bias would be evaluated by the sampling 
strategy and, subsequently, dumped. However, if the spacecraft 
were not in solar inertial before , which is the case 

for the Z-LV pass, then the solar inertial sampling strategy 
is not applicable. For this case, the momentum vector is 
sampled at Qp and the amount that is dimmed from (5) 
to (5) is the amount by which at (?) differs from 

the desired bias value [4] j. Note that, if this dump is 
executed exactly, at (?) is equal to the desired value. 

In a similar way, when thespacecraft is in solar inertial 
for the orbit preceedmg (5) , the level of (0) which 

determines the value of at , can be established with 

the existing solar inertial sampling strategy. The details 
of this procedure coinprise the contents of Appendix A. Thus, 
the importance of (?) is that, since it can be arranged to 
occur at the end of a dump cycled, its value can be specified 
by existing ATMDC software. For this reason it is the starting 
point for this study as it was for [1] . 

Orbital position (l) locates the point where the 
maneuver to acquire the Z-LV attitude is begun [1] . For the 
1 g 1 <.30° constraint used in [ 1] , point always occurs 

during the dark portion of the orbit. Point © must be 
selected some time after effective spacecraft sunrise to 
allow: correction of any attitude errors that may have 

accumulated on the dark side, and the updating of the strap- 
down calculations before beginning the maneuver. The angle y^ 

in Figure 1 was made variable because the duration of the dark 
portion of the orbit varies with g. For this study, y 

was taken to be y-^g) and point © taken fixed with respect to 

effective sunrise. In [1] , y^ was taken to be constant at the 

value for g=0°. This matter is discussed m more detail in 
§111 of this report. 


*This is also the biasing that was used in [1] . 

tThe regular solar inertial dump cycle is performed over an 
entire dark side pass rather than over 90 orbital degrees. 
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The rest of the sequence proceeds in a straight- 
forward manner (See Figure 1) . At , the acquisition rate 

is terminated and the orbital rate is initiated. Between 
and the spacecraft is in the Z-LV attitude. At ( 3 ) , tKe 

orbital rate is terminated and the acquisition rate for the 
return (to solar inertial) maneuver is initiated. When this 
rate is terminated at ( 4 ) , the spacecraft has reac qu ired 
solar inertial. It is held in this attitude until , 

when H is measured and the gravity-gradient dump maneuver 
begins . 


In [1] , the Z-LV maneuvers are svmmetrical 
about orbital noon (0=0). That is: and are equidistant 

from 6 = 0 ; = Y 4 / an ^ = Y 5 . The results of a 

digital computer simulation of such a symmetrical sequence 
combined with the biasing of defined by ( 2 ) are shown in 

Figure 4. In addition, two other runs were made. These are 
illustrated m Figures 5 and 6 and differ from the run shown 
m Figure 4 in that the solar incidence angle, 3, is taken to 
be 30° and 50°, respectively.* Note that the component of CMG 
angular momentum along the Y-geometnc axis, not the axis 
normal to the orbital plane, is the one that is zeroed in 
the latter two runs. This is because in an orbit with a Z-LV 
pass, the variation of (subscript Y denoting the geometric 

Y axis) exceeds in magnitude the variations of H and H , and 

A Z 

because the existing ATMDC software can only bias H relative 
to geometric axes. Also, as discussed in §ill, the level of 
projection of H onto the geometric X and Z axes, i.e. and 

H z , are also set to zero at as they were in [ 1] . The 

TACS impulse expenditure for these three cases is 20 lb-sec, 

28 lb-sec and 24 lb-sec, respectively.^ These results do not 
differ significantly from the results of [ 1 ] . 

III. An Operational Procedure for Biasing H y for Symmetrical 
Passes 


The momentum profiles m Figures 4 through 6, combined 
with the existing capability for biasing the level of H y (as well 


*The 8=50° run, while violating one of the ground rules of 
[1], is included because of the aforementioned change in the 
baseline that occurred after the publication of [1] . 

+ The minimum impulse bit used for this study is 4 
lb-sec. This number was inferred from section 13 of [4] . 
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as and H ^) , suggest an approach for reducing the amount of 

TACS propellant consumption even further. No£e from these 
figures that large changes m Hy occur from to © • 

Further, as explained in [6] , a substantial component of bias 
momentum is built up about the axis normal to the orbital 
plane when the spacecraft is m an Earth-pointing attitude. 

The sum of these two effects produces the largest change in 
H during the orbit. Thus, if the total change in H from 
to ( 3 ) can be centered between the TACS firing thresholds 
estaElished by the nested control law, then the entire EREP 
sequence can be executed by the CMGs with no TACS expenditure. 
This change from © to (3j is calculated by adding the 
following changes. 

1. AH due to termination of maneuver rate. 


AH Y21 = *£y ft- lb-sec, (4) 

2 

where: I = inertia matrix, slug-ft , 

= acquisition rate, (sec) for the eigenaxis 
maneuver from solar inertial to Z-LV,* 

a y = unit vector along the positive Y-axis of the 
spacecraft. 

2 . AH due to initiation of orbital rate. ' 

iH Y22 = “ ( I w o ) * ft- lb-sec, 
where: uj 0 = orbital angular rate, (sec) 1 . 

3 . AHy due to bias momentum accumulation in Z-LV attitude. 

2 

AHy23 = 3(0 o I xz At f t- Ik-sec. 

where I xz = the 1, 3 element of the inertia matrix in 

slug-ft 2 , and 

At = the total time spent in Z-LV attitude. 

For 0=0 0 and an Experiment Pass Region (EPR) defined 
by y 3 = 120° (See Figure 1) , this total change exceeds the 

range permitted by the nested strategy for no TACS firing. 


*In spacecraft coordinates, = 
positive about the -Y geometric axis. 


[0 -<^ 0] . That is wj is 
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However, some propellant can still be saved by positioning 
in Figures 4 through 6 as closely as possible to the TACS firing 
threshold.* Therefore, whether or not any TACS firing will be 
required, it is possible by simple calculation to obtain a 
desired value for H just before the termination of the EREP 

acquisition maneuver rate. This value, m turn, corresponds 
to a unique value of H at (5) , which can be established m 

three steps. 

1. Integrate the Y component of H in equation (5) m 
reverse time from © to © • 

H + u 1 X H as T^ m - X (5) 

where: is the vector representation, m spacecraft co- 

ordinates, of the maneuver rate and T^ m is the analytical 

representation of the gravity-gradient torque during the 
acquisition maneuver. 

2. Calculate the change in required to initiate the 

acquisition rate and subtract it from the final value of the 
previous integration. That is, if (t-^+) defines the 

desired value of momentum immediately after acquisition of 
the maneuver rate, and H y (t^-) is the value of momentum 

immediately before this rate acquisition, then 


H Y ^1 ^ H Y ^ t l + ^ + AH Y21 
where AH^^is defined by (4) . 

3. Integrate (6) in reverse time from ® to © . 

K - T ys (6) 

where T Yg is the Y component of the gravity-gradient torque in 

the solar inertial attitude. The value of H„ calculated by the 

Y J 


*The threshold used for this study was 6412.5 ft- 

Ib-sec. This was calculated using a conservative value of 2250 

ft- lb-sec for the individual CMG angular momentum (with the new 

wheel speed) in the equation, H =(.95) (3) (2250) = 6412.5 ft- lb- 

max 
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latter integration is the desired value of Hy at . As 

previously mentioned, a description of how this value can 
be achieved with existing ATMDC software is explained in 
Appendix A. 

Note that there is no need to perform the above 
calculations on-board. The desired value of Hy at Qy could 

be calculated on the ground and up- linked 4 or 5 orbits* 
before the scheduled Z-LV acquisition orbit. 

The above procedure cannot be extended in a simple 
way to establish optimal values for H and H . This is 

A Is 

because these axes are rotating relative to an inertial co- 
ordinate system** during various segments of the orbit 
described m Figure 1. Thus, no single value of H^. and 

can be chosen to minimize the magnitudes of their variation. 
However, very little performance penalty accrues from not 
biasing H and H„ since their variations are much smaller than 

the variation of Hy. For the computer runs performed in this 

study, H^. and were set to zero at © and produced good 

results . 

The result of applying this procedure to a case 
that is similar to the one shown in Figure 4 i.e. a 120° 
symmetrical experiment pass region and $ = 0°, is Hy (8 0 ) = 

390 ft-lb-sec. The momentum profile that obtains is plotted 
in Figure 7. The TACS impulse expenditure for this case is 
8 lb-sec as opposed to the 20 lb-sec expended if Hy(e c ) =0. 

Figures 8 and 9 show the resulting momentum profile for the 
same and Hy (0 O ) = 390 ft-lb-sec, except that 8 = 30° 

and 50°, respectively, for these two cases. The corresponding 
TACS impulse expenditures are 8 lb-sec and 4 lb-sec correspond- 
ing respectively to 28 lb-sec and 24 lb-sec expended when 
Hy (0 O ) =0 (Figures 5 and 6) . Note that even better results 

could have been obtained for 8^0 cases if a new evaluation 
of the reverse integration had been made rather than using 
the desired value for 8 =0 case. 


* This number, a result of simulations run at MSFC, was 
obtained on a telephone conversation with H. F. Kennel on 
November 30, 1970. 

**An exposition of the variation of H in a rotating 
spacecraft is given m [6] . 
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Actually, the plots in Figures 7 through 9 differ 
slightly from their counterparts in Figures 4 through 6 in 
two additional ways. Namely, they use variable values of 
Y ^ and y^ in Figure 1. These modifications, discussed below, 

could be added to any ground-based computation of the desired 
value of at . 

The variability of y-j_ is possible because of the 

variation of the duration of the dark portion of the orbit 
with £. If the dark portion of the orbit is greater than 
90 orbital degrees, which insures completion of a dump 
maneuver, then (T) can be selected with respect to effective 
spacecraft sunrise. This allows a longer maneuver time for 
cases when |ej > 0. For this study, (T) occurs approximately 
one orbital degree after effective sumfise. 

The variability of Y 4 can be used to provide an 

additional CMG momentum dump capability about the Y-geometric 
axis. That this capability may be of some importance is seen 
by noting that the value of at (T) in Figure 4 is approxi- 
mately -1500 ft- lb-sec. Moreover, this quantity was calculated 
by considering only the effects of the gravity-gradient torque. 

If the aerodynamic and venting torques were included, (the 
latter having a relatively large effect on the Y-axis) H y 

could be increased substantially. Depending on the size of 
the increase, could achieve a value that could not be dumped 

simultaneously with the X and Z axis momentum dump requirements. 
The details of a procedure for providing extra Y-axis dump 
capability are given m Appendix B. 

IV. Biasing for Unsymmetrical Passes 

As mentioned in §11, the nested control system must 
have the capability to permit performance of the EREP experi- 
ments during some 60° segment of the symmetrical region defined 
by y 3 = 120° in Figure 1. Thus, it is possible to produce passes 

that are not symmetrical with respect to orbital noon. The most 
unsymmetrical passes are: 

1. when © is 60° before orbital noon and ( 3 ) is coincident 
with orbital noon. 

2. when ( 2 ) is coincident with orbital noon and ©. occurs 
60° later. The first of these two cases will be called the noon 
exit case and the latter the noon acquisition case. 
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As shown in Figure 10, the noon exit case can be 
handled by the same procedure used for symmetrical passes 
with an added 8 lb— sec/orbit impulse penalty (16 lb— sec 
total) . If H y is zeroed at (o) , the impulse expenditure 

is 24 lb-sec. Either of these may be considered acceptable. 
However, the setting of H y to zero at Q0) for the noon 

acquisition, 6=0° case results m a 44 lb— sec TACS impulse 
expenditure (see Figure 11) . Further, not much improvement 
occurs if H y at @ = 390 ft- lb-sec. Thus, some generalization 

of the procedure is necessary if unsymmetrical passes are to 
be considered. This is done as follows 


1. Determine AH y (t 2 ) = AH y21 +AH Y22 at(^as in §111. 

2. Determine AH y (t 3 ) = AH y at as follows: 

a. AH y due to termination of orbital rate 

AH y3 i — { I_to <> ) • a y ft— lb— sec — — AH y 22 * 

b. AH y due to initiation of return maneuver rate 

AH Y32 = ^—2^ * — y 


where ^ — the vector representation in spacecraft co- 
ordinates of the return maneuver rate. 


3. If J AH y (t^) | <J AH y (t 2 ) | , use the symmetrical procedure 

in Sill. That is, calculate AH y from ( 2 ) to ( 3 ) and center, 

integrate (5) in reverse time, calculate AH y2 ^,and integrate (6) 
m reverse time. 


4. If 


AH y (t 3 ) 


> AH y (t 2 ) 


proceed to 5. 


5. Center AH y (t 3 ) with respect to the thresholds 

established by the nested control system strategy. The 
greatest lower bound of the centered profile provides the 
boundary value for the following integration. 


6. Perform a reverse time integration on the equation. 
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v 


where T „ = the gravity-gradient torque acting on the Earth - 

— gi, 

pointing spacecraft. This integration is carried out from 
to © ( see e.g., Figure 11). 

7. Algebraically add the negative of AHy(t 2 ) to the 
final value of the integration performed m 6 . 

8. Continue as with the symmetrical procedure, i.e. 

integrate in reverse time from © to , calculate AH^Ct^) 

and add it appropriately to the final value of the previous 
integration, and then integrate m reverse time from Q) to (o). 

When the noon acquisition, 3 = 0° case was run with 
the value of H at specified by the above calculation 

(i.e., -2610 ft- lb-sec) t the TACS impulse expenditure was 
reduced from 44 lb-sec to zero lb-sec (see Figure 12) . 

V. Conclusions 


While the results of [1] are adequate within the 
framework of their groundrules, i.e., symmetrical Z-LV passes and 
|3[<30°, the possibility of a more ambitious EREP warrants 
the development of an operational technique for biasing H. 

The approach taken m this memorandum was to extend the 
results of [1] without trying to change the baselined 
maneuver sequence.* Three extensions to the results in [1] 
were obtained: 

1. An operational procedure was given for determining 
a value of ( 0 0 ) that reduces the TACS expenditure for 

symmetrical Z-LV passes. 

2. An extension of the above procedure for unsymmetncal 
maneuvers was developed. 

3. Some observations were made in Appendix B on providing 
extra dump capability for the purpose of extending the number 
of successive Z-LV passes that can be made without increasing 
TACS expenditure. 

Finally, it should be mentioned that the above 
operational procedure for determining H Y (0 O ) is of limited 

value if only two CMGs are operative. In such a case recourse 
might be made to the maneuver sequences described in [6] for 
substantial reduction of TACS expenditure. 


*An analysis of other possible maneuver sequences was 
presented m [6] . 

tThis depends on the ability of the CMGs to hold the space- 
craft in solar inertial with, say, a 1305 ft-lb-sec bias and no 
TACS firings. This seems reasonable. 
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APPENDIX A 


SOME DETAILS ON THE SAMPLING TECHNIQUES OF THE 
CMG MOMENTUM DUMP SCHEME 


The reason for initializing the H vector at orbital 
position is imbedded in the details of how the CMG momentum 

is sampled at various points in the orbit. This sampling is 
performed for the purpose of determining how much CMG momentum 
should be dumped during the succeeding night-time pass. A 
complete description of the momentum dump scheme is given in 
Section 10 of [5] . The only aspects of that scheme that are 
presented here are those necessary for a description of how the 
ATMDC software would command the H vector to be positioned 
properly at . 

A general mathematical representation of the CMG 
momentum variation along the body axes is 


H i (0 ) = Hq^ + (0+a^— 3ir/4 ) + H 2 ^cos 2 (0+a^) 


(A-l) 


where: 0 = w^t = the orbital angular excursion measured from 

orbital noon, = the aforementioned angular offset due to the 

non-parallel relationship between principal and geometric axes , 
i = x,y,z, H . = constants and j = 0,1,2. Note that the 

momentum variation is comprised of a constant term, a term that 
increases linearly with time or orbital angle relative to , 

and a sinusoidally varying term. The constant term, H^, estab- 
lishes the CMG momentum level and is the quantity that can be 
adjusted by the dump regime. As previously mentioned, the 
following definitions can be made: 

1. dump midnight, (6) , is defined by 


= -mir , m = 1,3,5,... 


0 + a 


r 
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2 . dump noon is defined by 0 + a = 0 ,* 

3. H, (n) is the value of H measured at a point kiT/4 

from dump midnight, k = 2,3,4, 5,6, 7, during the n — 
consecutive orbit containing a Z-LV pass, n - 1,2,... . 


To complete the preliminaries , the assumption will be made that 
— 0 m (A-l) . This assumption is nominally true because of 

the Z-axis rotation that is selected to minimize the bias 
momentum accumulation about the Y and Z axes of the spacecraft. 
The Y-axis momentum variation will, then, be described by 


H ( 0 +a ) = H q + H 2 cos2(9+cc) 


(A- 2) 


Note that the Y subscripts have been dropped. 

The ATMDC equations for the amount of momentum to be 
dumped to achieve a certain level for H( 6 ) is given by (A-3). 

In this equation 

H A (n) = |[H 3 (n) + Hj. (n) ] - H B (A-3) 

\ 

H Q (n+l) = H Q (n) - H A (n) . (A-4) 

Hg is an input that corresponds to the desired level of the 

momentum variation and to the desired value of H^. at . 

The connection between the level and H^fe^) will be shown m 

the next paragraph. Equation (A-4) expresses the fact that the 
level of H (i.e., H n ) during the n + first orbit equals the 

u th 

level of during the n — orbit less the amount dumped during 

*The = 0 ^( 6 ) are slightly different for each axis. We use 
here and from now on a — since the initialization procedure is 
not performed for the X and Z axes . 
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the n — damp cycle. This assumes that the actual dumped momentum 
is equal to the commanded dump momentum. 

The sampling instants m (A-3) , i.e., k - 3,5, have 
been chosen to remove the effect of the sinusoidally varying 
term in (A-2). To show this, note that. 


H 3 (n) = H Q (n) + H 2 (n)cos2 (-■ ir/4) = H Q (n) 


H 5 (n) « H Q (n) + H 2 (n)cos2 (tt/ 4) = H Q (n) 
which when substituted into (A-3) yields 

H (n) = H (n) - H B . (A-5) 


(A-5) yields the amount of H Y to be dumped for the purpose of 

establishing a new level*. The value of this new level is 
established by (A-4) , 


Hq (n+1) = H 0 (n) - (n) = H Q (n) - H Q (n) + Hg = H fi . (A-6) 


Note that this level will be maintained until a desired value 
other than Hg is inputted. This is shown by solving (A-3) and 

(A-4) for H^(n+1) and Hq (n+2) . 


H A (n+l) = H 0 (n+1) - H B = 0 


Hq (n+2 ) = Hq (n+1) - H A (n+l) = H £ 


*This momentum has been separated from the momentum that 
will be dumped, e.g., about the X-axis, due to such terms in 
(A-l) as H lx ? 0. 
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Finally, it is easy to show that H y at of the n + first 

orbit after n orbits m solar inertial is equal to Hg. To do 

this, note that © corresponds to (6+ct) = - 3ir/4 . Thus, 

(A- 2) becomes 


H (n+1) = H q (n+1) + H 2 cos (-3tt/2) , 


If / 0 in (A-l) , a digital filter equation [5] 

governs the convergence of to its desired level. This level 

will not, in general, be reached after a single dump cycle. 
Therefore, as mentioned in §111, the off-line calculation of 
Hg should be up-linked 4-5 orbits before the schedule Z-LV pass 

to allow enough time for the system to readjust the momentum 
level. 
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AN ADDITIONAL Y-AXIS GRAVITY- GRADIENT DUMP CAPABILITY 


As mentioned in §111, an extension of the presently 
baselined maximum of two successive Z-LV passes raises the 
question of whether or not the momentum dump cycle can achieve 
successively new values of at . Were this to occur, the 

desired value of at © could not be attained and more than 

the minimum amount of TACS propellant would be consumed. The 
purpose of this Appendix is to show how to deal with this situ- 
ation by varying y^ (see Figure 1) to provide an additional 

Y-axis gravity-gradient dump capability. 

B-l A General Description of the Procedure 

Consider once again, the momentum profile in Figure 4. 

Note that, for a symmetrical pass, the momentum variation would 
be symmetrical were it not for the bias term and the TACS firings 
that occur during the Z-LV attitude hold segment. This asymmetry 
can be (partially or totally) compensated for by making T 4 < ^2 

(see Figure 1) . A heuristic explanation of this statement is 
based on two principles : 

PI. the sinusoidal variation of the gravity-gradient torque 
as shown in Figure 2 , and 

P2. the inverse relationship between the momentum variation 
due to the gravity-gradient torque and the rate of rotation relative 
to a set of local vertical coordinates . * 

It is also helpful to note that, during both the acquisition and 
the return maneuvers, the spacecraft moves at a rate w > rela- 
tive to local vertical coordinates. 

Refer to Figure 1. From © to the spacecraft is 

in solar inertial and, for the case shown, is always in a region 
of negative gravity-gradient torque about the Y-axis. Similarly, 
from to © , the spacecraft is in solar inertial and, for 

the case shown, always in a region of positive gravity- gradient 
torque about the Y-axis. If y^ = y,-, the net change in over 


*An example of P2 is given in §B-2 of this Appendix. 
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these two segments is zero. In the same fashion, a look at 
Figure 4 shows that there is a net positive change from 
to © (neglecting that AH^ for maneuver acquisition) and a 

corresponding negative change from to ^ (also neglecting 

AH^. for maneuvers). Since the bias accumulation is negative, 

from ( 2 } to © , it would be advantageous if the net positive 

AHy on the afternoon side of the orbit exceeded the net negative 

change on the forenoon side. That this can be accomplished by 
setting y 4 < y 2 ( an d therefore, y 5 > Y^) is a direct result of 

and P 2 - The mathematical statement of this procedure is 


6H 01 + 6H 12 + 6H 34 + 6H 45 > 0 ' 


where 6 H. . = H. - H . Whether or not this increase in H„ is 
13 3 1 Y 

enough to offset the negative accumulation during the Z-LV pass 

depends on the existence of a solution for y 4 in (B-l) that 

does not violate momentum constraints. 


6 H n . + 6 H no + 6 H 0 _ + 6 H_ + 6 H-. + 6 H. C = 0 

01 12 23 f 34 45 


(B-l) 


where = the amount of momentum put into (or taken out of) 

the system by TACS firings. 

There are three points to be noted about (B-l) 

1. It does not take into account the AH^'s necessary 

for the acquisition and removal of various maneuver rates. This 
is because, if no TACS assist is used for these functions, they 
cancel each other. In a case such as the one shown in Figure 11, 
(B-l) would have to be modified to include these AE^'s. 

2. The determination of y 4 m (B-l) represents a 

solution to a transcendental equation and must, therefore, be 
done numerically. An example of a calculation of the solution 
to (B-l) for a particular case is given in §B-2 of this Appendix. 
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3. If the solution of (B-l) implies a rate that 
violates momentum constraints, then one of two options could 
be used. 

a. The y 4 that corresponds to the largest allowable AH^ 

for maneuver acquisition could be used. This was done for the 
run illustrated by Figure 7 and produces an effective dump of 
600 ft-lb-sec. Note the improvement relative to Figure 4. 

b. Another angle, say y 3 , could also be made variable. 

While this makes (B-l) a single equation with two unknowns, it 
might be possible to incorporate a constraint equation, e.g., 
the momentum constraint, for the purpose of achieving a unique 
solution. 

B-2 An Example 

In §B-1, a description of a procedure for providing 
additional Y-axis dump capability was presented. The calculation 
to be made was expressed by a functional equation (B-l) , which is 
repeated below for convenience. The purpose of this section is 
to provide an example of the type of calculation required to 
solve (B-l) . 


6H q1 + <SH 12 + 6H 23 + 6H f + SH 34 + 6H 45 = 0 (B-l) 


where = HL - = H ( j ) - H( i ) , i = 0 ,1 ,2, 3, 4; j = i+1. 

The solution of (B-l) will be calculated for the run 
that is illustrated in Figure 7. Since y 4 is the only variable, 

the first four terms of (B-l) are constant. The functional forms 
corresponding to these terms are easily determined and will be 
omitted. From Figure 7, 



6H 01 

H 1 " 

H 0 

= -2321 

- 390 = -2711 



6H 12 = 

H 2 - 

H 1 

= 6399 

- 3744 = 2655 


6H f + 

6H 23 “ 

H 3 - 

H 2 

= -6408 

+ 4806 = -1602 


'• 6H 34 + 

6H 45 = 

1658 

- ' 

1 

H 

O 

HI 

6H 12 " 6H 23 " 6H f 

(B-2) 
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The functional form of <$H 34 1S determined as follows. 

If t^ and t^ are the respective times corresponding to © and 

® , then, for a spacecraft maneuvering at a rate relative 
to solar inertial , 


<SH 


34 


H, 



T (t)dt 
gr 


(B-3) 


where T = the Y-axis component of the gravity- gradient torque. 
This can be shown to be equivalent to 


6H 


7 085 

34 = 2(a) 0 +a) 2 T {cos2[{w O +(0 2 ) (t 4" t 3 )+al “ COs2a} ' 


(B- 


The unknown quantities in (B-4) are t 4 and However, these 

quantities are related by 


“2 (t 4 ' t 3 ) = tlJ 0 t 3 ' (B ~ 5) 

Substituting (B-5) into (B-4) , letting 6^ = o>gt 4 and 0^ = Wgt^, 
'and using Wg = 1.11845 x 10 ^ (sec) \ finally leads to 


4 ® 3^ 

6H 34 = 3167.3 - {cos2 (0 4 + a) - cos2a} 


(B— 6 ) 


where 0 4 is the only unknown. 

From to the spacecraft is in the solar inertial 

attitude. The gravity-gradient torque for this mode is illustrated 
in Figure 2 and is given by 


T y = -7.085 sin2 (oigt + a). 
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Therefore , 

7 n ft 'i * 

6H 45 ” {cos2 (bigt,. + a) - cos2 (to Q t^ + a)} . (B-7) 

For this example, cu^t,. + a = 3ir/4. This implies that 

6H^ = —3167.3 cos2 (0^ + a) . {B-8) 


The substitution of (B-6) and (B-8) into (B-2) yields 
a transcendental equation that must be solved numerically. For 
this example, the solution of (B-2) violated the momentum con- 
straints on the system imposed by the nested control strategy. 
The alternatives for this case were discussed in §B-1. 


*The terms outside the curly brackets m (B-4) 
illustrate P2 in §V. 


and (B-7) 



ORBITAL 

NOON 



ATTITUDE PROFILE 


FROM: 

® TO (O) , SPACECRAFT EXECUTES DUMP MANEUVER. 

(o) TO (T) , SPACECRAFT HELD IN SOLAR INERTIAL ATTITUDE. 

® TO (2), SPACECRAFT MANEUVERS TO Z-LV ATTITUDE. 

® TO ®, SPACECRAFT HELD IN Z-LV ATTITUDE 

® TO (4), SPACECRAFT MANEUVERS BACK TO SOLAR INERTIAL ATTITUDE. 

(4) TO @ , SPACECRAFT HELD IN SOLAR INERTIAL ATTITUDE 

(5) TO (0) , SPACECRAFT EXECUTES DUMP MANEUVER. 


FIGURE 1 - PROCEDURE FOR ACQUISITION AND HOLD OF Z-LV ATTITUDE 
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Ty - component of gravity-gradient torque acting about the Y geometric axis when it is 
parallel to the normal to the orbital plane. 

6 - orbital excursion angle, 6=0 corresponds to orbital noon, 6 = ±7T corresponds to 
orbital midnight. 

FIGURE 2 - GRAVITY-GRADIENT TORQUE VS ORBITAL EXCURSION 



Hy - MOMENTUM VARIATION PRODUCED BY Ty IN FIG. 2 


FIGURE 3 - CMG ANGULAR MOMENTUM VS ORBITAL EXCURSION 


FIGURE 4 - Y-AXIS COMPONENT OF H VS ORBITAL ANGULAR EXCURSION, 6 . 

120° SYMMETRICAL Z-LV PASS, j3 = 0° 

H' {(o)} =[000] = H'(0 O ), 2250 FT-LB-SEC PER CMG 
SPACECRAFT IN SOLAR INERTIAL FROM @ TO (O) 








(FT- LB-SEC) 




FIGURE 5 - Y-AXIS COMPONENT OF H VS ORBITAL ANGULAR EXCURSION, 0 . 

120° SYMMETRICAL Z-LV PASS, |3 = 30° 

Jf (©) =[0 0 0] = H'(0 o ),225O FT-LB-SEC PER CMG 
SPACECRAFT IN SOLAR INERTIAL FROMl© TO © 



FIGURE 6- Y-AXIS COMPONENT OF H VS ORBITAL ANGULAR EXCURSION 

120° SYMMETRICAL Z-LV PASS,/3 = 50° 

H' ((0)) = [ 0 0 0 ] = H' «? 0 ), 2250 FT-LB-SEC PER CMG 
SPACECRAFT IN SOLAR INERTIAL FROM ® TO ( 5 ) 





(FT-LB-SEC) 




FIGURE 7 - Y-AXIS COMPONENT OF H VS ORBITAL ANGULAR EXCURSION, 0 . 
120° SYMMETRICAL Z-LV PASS, ) 3 = 0° 

H' (®) = [0 0 0] = H'(0 o h 2250 FT-LB-SEC PER CMG 
SPACECRAFT IN SOLAR INERTIAL FROM (S)TO 1® 



(FT-LB-SEC) 




FIGURE 8 - Y-AXIS COMPONENT OF H VS ORBITAL ANGULAR EXCURSION, 6 . 
120° SYMMETRICAL Z-LV PASS, jS= 30° 

H'((o)) = I 0 0 0 ] = H' {0 o }/225O FT-LB-SEC PER CMG 
SPACECRAFT IN SOLAR INERTIAL FROM (§) TO (o) 



(FT-LB-SEC) 




FIGURE 9 - Y-AXIS COMPONENT OF H VS ORBITAL ANGULAR EXCURSION, 6 

120° SYMMETRICAL Z-LV PASS,0 = 50° 

H' ((0)) =[0 0 0]= H'(0 O >. 2250 FT-LB-SEC PER CMG 
SPACECRAFT IN SOLAR INERTIAL FROM'(5) TO ( 5 ) 


(FT-LB-SEC) 



FIGURE 10 - Y-AXIS COMPONENT OF H VS ORBITAL ANGULAR EXCURSION, 6 
120° SYMMETRICAL Z-LV PASS, 0 = 0° 

H' ((o)) = [0 0 0] = H'{0 o ),-225O FT-LB-SEC PER CMG 
SPACECRAFT IN SOLAR INERTIAL FROMj© TO ( 5 ) 





FIGURE 1 1 - Y-AXIS COMPONENT OF H VS ORBITAL ANGULAR EXCURSION, 6 
60° UNSYMMETRICAL Z-LV PASS (NOON ACQUISITION) 

H'((0)) = [0 0 0] =H'(0 o ),225OlFT-LB-SECPERCMG 
SPACECRAFT IN SOLAR INERTIAL FROM® TO ® 





(FT-LB-SEC) 
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FIGURE 12 - Y-AXIS COMPONENT OF H VS ORBITAL ANGULAR EXCURSION, Q. 

60° UNSYMMETRICAL Z-LV PASS (NOON ACQUISITION) 

H'((o)) = [0 0 0 3 = H'(0 O ), 2250 FT-LB-SEC PER CMG 
SPACECRAFT IN SOLAR INERTIAL FROM (?) TO (o) 





